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ABSTRACT 

We present the SLoWPoKES-II catalog of low-mass visual binaries identified from the Sloan Digital Sky 
Survey by matching photometric distances. The candidate pairs are vetted by comparing the stellar density at 
their respective Galactic positions to Monte Carlo realizations of a simulated Milky Way. In this way, we are 
able to identify large numbers of bona fide wide binaries without the need of proper motions. 105,537 visual 
binaries with angular separations of ~l-20", are identified, each with a probability of chance alignment of 
<5%. This is the largest catalog of bona fide wide binaries to date, and it contains a diversity of systems— 
in mass, mass ratios, binary separations, metallicity, and evolutionary states—that should facilitate follow-up 
studies to characterize the properties of M dwarfs and white dwarfs. There is a subtle but definitive suggestion 
of multiple populations in the physical separation distribution, supporting earlier findings. We suggest that wide 
binaries are comprised of multiple populations, most likely representing different formation modes. There 
are 141 M7 or later wide binary candidates, representing a 7-fold increase in the number currently known. 
These binaries are too wide to have been formed via the ejection mechanism. Einally, we find that ^6% of 
spectroscopically confirmed M dwarfs are not included in the SDSS STAR catalog; they are misclassified as 
extended sources due to the presence of a nearby or partially resolved companion. The SLoWPoKES-II catalog 
is publicly available to the entire community on the world wide web via the Eiltergraph data visualization portal. 
Subject headings: (stars;) binaries: visual — stars: low mass — stars: brown dwarfs — stars: late-type — 
(stars;) white dwarfs — (stars:) subdwarfs stars; statistics — 


1. INTRODUCTION 

Components of binary (or multiple) systems are ideal co¬ 
eval laboratories to study star formation, to benchmark stellar 
evolutionary models, and to calibrate empirical relations that 
determine fundamental stellar parameters. While detailed and 
precise measurements of individual objects or small samples 


provide important tests for evolutionary models (e.g., |White 
et al.||199^ Stassun et al.||2007 2008| l, large statistical sam¬ 


ples are necessary for properly constraining the behavior and 
intrinsic variation of star formation and of stellar properties. 
All of these rest on the premise that individual components 


in a stellar system are formed at the sa me time (White & 
|Ghez 2001[|IGaus & Hillenbrand|2009a|l, of the same mate 


rial ( [Schuler et al.|2011[|Dhital et al.|2012| l, and have evolved 
in the same environment. 

Until recently, the intrinsic faintness of low-mass stars— 
generally defined as the regime brack eted by the hydrogen 
burning minimum mass, ~0.075 Mq (jBurrows et al.||1997|l, 
and the onset of molecu lar lines in the photosphere, ~0.8 Mq 
(|K irkpatrick et al.|1991|l —^has limited studies to small, nearby 
samples. Almost two decades ago, the Palomar-Michigan 
State University (PMSU) survey cataloged the spectr a for 
^3,000 M dwarfs (dMs) and w as the largest such study ( |Reid| 
jet al.| 1993) [Hawley et al.|1996| l. Other samples used to study 
the distribution of low-rnass binaries and to calibrate low- 
mass stellar prop erties using binaries were even smaller, with 
30-100 systems ( Eischer & Marcy|1992[ Henry & McCarthy 
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|1993[[Reid & Gizis|1997[|Delfosse et al.|2004| l. Hence, these 
studies often lacked the statistical robustness for firm, inde¬ 
pendent results, which were then often tethered to and studied 
in comparison with higher-mass stars. In addition, model¬ 
ing efforts for the low-mass late-K and M dwarfs have been 
stymied by incomplete molecular line lists resulting in un¬ 
certain opacities caused by the molecules in the stellar pho- 


tosphere at effective temperatures of <4300 K (Hauschildt, 
[Allard, & Baronl|l999| l. As a result of these observational 
and modeling challenges, low-mass stars have not been well 
characterized, and our techniques for measuring their prop¬ 
erties are ill-defined. However, lo w-mass stars comprise of 
>70% of the stars in the Galaxy (Henry 1998 i and are the 
best tracers of its distributio n and (at least nearby) structure 
(e.g. [Bochanski et al.pOlOj l. Low-ma ss stars also have life ¬ 
times longer than that of the Galaxy ( Laughlin et al.|[T997 l, 
making them the ideal tracers of its formation, chemical, and 
dynamical history. 

The advent of deep, all-sky surveys has revolutionized low- 
mass star and brown dwar f science. In par ticular, the Sloan 
Digital Sky Survey (SDSS;|York et al.|2000]l, the Two Micron 
All Sky Survey (2MASS; |Cutri et al. 1^0031), and UKIR T In¬ 
frared Deep Sky Survey (IJKJDSS; Lawence et al. 2007 1 , and 
Wide-held Infrared Survey Explorer (e.g., [Wright et al.|2010 1 
have been critical in expanding sample sizes. SPSS alone 


has enabled a photometric catalog of >30 million ( [BochansE 
et al. 20101 and spectroscopic catalog of >70,000 |West et al. 
2011 hereafter, Wll) low-mass stars. [Munn et al.| ( [2004'[ 
2008 I combined the USNO-B and SDSS astrometry to cal- 
culate proper motions for SDSS photometric catalogs. This 
proper motion catalog is 90% complete down to p < 19.7, 
with typical errors of 5-7mas yr“^. However, the lower reso- 
lu tion limits the catalo g to sourc es separ ated by ^ 7". 

Dhital et al.[ (120101 hereafter, [Paper Ijl identihed the Sloan 
Low-mass Wide Pairs of Kinematically Equivalent Stars 
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Figure 1. The average number of visual companions (left) and real binaries (right), as ass essed by our Galactic models, as a function of angular separation in 
>1600 lines-of-sight where binary candidates were identified in the SLoWPoKES catalog jPaper IJ. The solid histograms shows the SDSS distributions while 
the dashed histograms show the distributions generated in our Monte Carlo models. In this figure, we show the distribution of the average number of companions 
around >1600 candidates. We succesively match the angular separation (with the results shown in red histograms), photometric distances (blue), and proper 
motions (purple). The model distributions predict more stars at larger angular separations as they include all stars whereas the data count only the stars that are 
detected in SDSS. Intriguingly, the blue and purple histograms were almost identical up to 20-25”, suggesting proper motions were not needed to identify those 
binaries. This is seen even more clearly in the right panel. In this paper, we conduct a search for binary companions around >^24 million low-mass stellar sources 
in the SDSS survey without using proper motions. 


(SLoWPoKES) catalog of wide binaries by matching posi¬ 
tions, photometric distances, and vector proper motions. The 
fidelity of each parr was assessed by a six-dimensional Galac¬ 
tic m odel, which was built using empirical stell ar number den¬ 
sity ( Juric et ar]|2008 Bochanski et al, [2010 1 and space ve- 
locity distributions (|Bochanski et al.|2007|l. With 1342 com¬ 
mon proper motion pairs, with <5% probability of chance 
alignments, SLoWPoKES was the largest catalog of low- 
mass, wide binaries and contained a diverse set of pairs with 
G/K+dM, white dwarf (WD)-bdM, dM+dM, and M subd¬ 
warf (sdM+sdM) systems. SLoWPoKES has enabled a vari¬ 
ety of follow-up obse rvations to p robe higher-order m ultiplic¬ 
ity (Law et al. 2010|l, metallicity (Dhital et al.|2012|l, activity 


( Gunning et al.lpolT i, rotation, and the age-activity relation 
(Morgan et al., in preparation, Massey et al., in preparation) 
of M dwarfs. In addition, ( [Andrews et al.|2012| |2015] l have 
adapted the technique used for SLoWPoKES to identify dou¬ 
ble the sample of known wide WD-fWD binaries. 

The SLoWPoKES catalog was restricted at small angular 
separations [9 < 1") and at lower masses (< M6) by its 


dependence on the SDSS/USNO-B proper motions (Munn 
jet al.|20d4) l. Specifically, due to the shallow faintness limit of 
USNO-B, a significant fraction of mid-late dMs that are de¬ 
tected in SDSS do not have USNO-B counterparts and, there¬ 
fore, proper motions. Similarly, the resolution of USNO-B 
sets a separation limit of ~ 7" for the SDSS/USNO-B proper 
motion catalog. These incompleteness are inherited by the 
SLoWPoKES catalog. In addition, using proper motions to 
identify binaries precludes systems that are either at large dis¬ 
tances or nearby but moving slowly with respect to the Local 
Standard of Rest. Eor example, we used a minimum proper 
motion of 40 mas yr“^ for the SLoWPoKES catalog. 

Results from the Galactic model indicated that visual bi¬ 
naries with small angular separations could be identified at 
a high level of fidelity matching photometric dis tances 
alone ([Paper Ijl. Figure [T] (left) shows Figure 5 from [Pap'CT 


|Ij The solid red histogram shows the number of stars, 
eraged among gtrsiml6QQ lines-of-sight where binary can¬ 
didates were identified, as a function of angular separation, 
as measured in the SDSS DR7 photometric catalog. The 
dashed red histogram shows the distribution for single, non- 


associated stars in the same fields as simulated by our Galac¬ 
tic model. In essence, this was a measure of the likelihood 
of whether the companion stars were real binaries or simply 
chance alignments. The blue and purple histograms show the 
distribution when the photometric distances and proper mo¬ 
tions were matched. Each addition of a dimension resulted 
in a rejection of a large number of chance optical pairs. The 
larger stellar number counts at large separations in our simula¬ 
tions (as compared to the data) was caused largely by the fact 
that the data were limited to r < 20.5 as the SDSS/USNO- 
B proper motions were required whereas our Galactic model 
simulates all stellar objects. Two features stood out in this 
figur e: (1) the exces s of pairs at small separations, which was 
how [Michell[ ([1767[l first identified binary systems, or “dou¬ 
ble” stars as he called them and (2) the blue and purple his¬ 
tograms were essentially the same until about 9 ^ 20 — 25", 
suggesting those binaries could have been identified without 
proper motions. This is even more evident in the right panel 
of Figure [T] where we have plotted the difference between the 
data and model distributions. These figures clearly demon¬ 
strate that proper motions are not required to identify bina¬ 
ries up to a critical separations with a high level of fidelity. A 
caveat is that the critical angular separation obviously depends 
on the magnitude limit of the sample. For the SLoWPoKES 
sample with a limiting magnitude of r = 20.5 and proper 
motion >40 mas yr“^, >90% of the binary candidates within 
20" with matching photometric distances also had matching 
proper motions and were classified as CPM pairs. 

In this paper, we extend the SLoWPoKES catalog by identi¬ 
fying binary systems with angular separations of 1-20" based 
entirely on SDSS photometry and astrometry. This allows us 
to identify visual binaries to r = 22.2, with significant num¬ 
bers at the mid-late M spectral types. In Sectionj^we describe 
the i nitial sam ple of low-mass stars that we search around. 
As in [Paper Ij our search algorithm is based on matching an¬ 
gular separation and photometric distances supplemented by 
a Monte Carlo-based Galactic model, which is described in 
Section [3] We discuss the characteristics of the resultant bi¬ 
nary sample in Section We examine the currently debated 
formation theories for wide stellar binaries and VLM/BDs in 
light of SLoWPoKES-II sample of binaries in Section]^ We 
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summarize our results in Section |6] 

The SLoWPoKES and SLoWPoKES-II catalogs, along 
with followup spectra, are publicly available onlineH 

2. SDSSDATA 

One of the largest and most influential astronomical surveys 
ever conducted to date , SPSS is a compr ehensive imaging and 
spectroscopic survey (|York et al.||2000|l. Over eight years of 
operation between 2000-2008, it collected imaging and spec¬ 
troscopic data for over a decade years using a dedicated 2.5-m 


telescope a t Apache Point Observatory, New Mexico ( Gunn 
|et al.||2006| . The telescope has a 120 meg apixel camera that 
has a field of view of 1.5 square degrees ( |Gunn et al.||1998| l 
and conducts imaging in five broad optical bands (u qrlz) be- 
tween ~3,000 and 10,000 A (Fukugita et al. 19961. The last 
data release with imaging data. Data Release 8 (DR8), com¬ 
prised of ^450 million unique objects over 14,555 square de¬ 
grees of the sky, spanni ng the entire norther n sky as well as the 
Southern Galactic Cap ( |Aihara et al.|201 l]l. The global abso¬ 
lute astrometric precision was 70 mas ( |Pier et al. 12003] ) while 
the photometry has relative calibra tion accuracies of 2% in th e 
u band and 1% in the griz bands ( jPadmanabhan et al.|2008| l. 
The catalog is 95% complete for point sources of 22.0, 22.2, 


et al 


22.2, 21.3, and 20.5 in the ugriz bands, respectively (Gunn 
.|19981l. The corresponding spectroscopic survey has ob- 


tained ~1.8 million optical spectra, with A/A A « 2000, over 
9274 square degrees (|Aihara et al.|2011|| using a pair of fiber- 
fed double spe ctrographs. ITie Third SToan Digital Sky Sur¬ 
vey (SDSS-III;|Eisenstein et al.|20111l, which is comprised of 
four different spectroscopic surveys, is currently underway. 

The DR8 photome tric catal og has more than 200 million 
point sources. As in [Paper T| we used the Catalog Archive 
Server query tool (CasJob^jito select the sample of low-mass 
stars from the DR8 STAR table as having r — i > 0.3 and 
i — z > 0.2, co nsistent with spectral types of K5 or later 
( West et al.|2008| l. Selecting from the STAR ensures that the 
object is a PRIMARY and not a duplicate detection and that 
its morphology is consistent with being a point source. This 
morphological class ification is >95% accurate up to r « 21.5 
( |Lupton et al.||200T|). Even up to r « 22.5, stars outnumber 
galaxies >8:1 (jPadely et al.|2012|l, and, therefore, we do not 
expect compact galaxies as frequent interlopers in our sam¬ 
ple. To ensure excellent photometry, we made cuts on the 
standard quality flags—all of PEAKCENTER, NOTCHECKED, 
PSE_ELUX_INTERP, INTERPXE NTER, BAD_COUNTS_ERROR, 
SATURATED were set to be 0 ( Bochanski et al.||2OT0t [Paper] 

—and required the errors in PSE magnitudes to be < 0.10. 
Throughout our analysis, these quality cuts were performed 
only for the bands that were used in the analysis for that par¬ 
ticular star: iz for M7 or later stars, riz for K5-M7 stars, griz 
for F0-K5 stars, and ugriz for white dwarfs. While no faint¬ 
ness limits were speciflcally adopted, the requirement on the 
error in PSE magnitudes effectively limits resulting samples 
to the SDSS 95% completeness limits. 

The resulting sample yielded 33,589,670 ste llar sources, 
with c olors consistent with K5-M9 dwarfs. In IDhital et al.l 
( |2010[ ) we found that chance alignments were unacceptably 
high near the Galactic Plane due to the higher stellar den¬ 
sity and the higher uncertainties caused by higher extinction. 
Therefore, we rejected stars at Galactic latitudes less than 20°. 
We also rejected stars with photometric distances larger than 

^ http://slowpokes.Vanderbilt.edu 

* http://skyserver.sdss3.org/CasJobs/ 


2500 pc as the distance uncertainties are larger than ^350 pc, 
making dista nce matching meaningless at those distances (see 
Section |3.2| for further discussion). As a result, we started 
with an initial sample of 24,036,982 low-mass stars., around 
which to search for companions. 

3. METHOD: IDENTIEYING BINARY CANDIDATES 
3.1. Assessing the SDSS source classification algorithm 

The ^hest binary we were able to identify was >1" (see 
Figure H below), suggesting a fundamental limitation in our 
technique at that value. However, with a plate scale of 
(y!396 pixel“^, SDSS should be able to distinguish tighter bi¬ 
naries. Both unresolved and partially-resolved binaries, with 
separations smaller than the resolution limit or the plate scale, 
have been ident ified from both the 2MASS ( [Kraus & Hil- 
lenbrand|2007b|l and the Palomar Transient Factory ([Terziev| 
et al.[2013|l photometric catalogs, albeit aided by multi-epoch 
data in the latter case. Therefore, we conducted an examina¬ 
tion of (1) what was restricting us from identifying binaries at 
<1" and (2) what fraction of resolved binaries we were miss¬ 
ing in our resultant sample due to this limitation. 

First, we briefly describe the SDSS point source classifica¬ 
tion scheme. The algorithm is fully documented on the SDSS 
webpage^and is beyond the scope of this paper, thus we only 
outline the scheme that classifies identified sources into the 
STAR table. The SDSS pipeline uses RESOLVE to identify 
unique sources and and DEBLEND to deblend them when 
multiple sources are present. All of the identified sources, for 
which photometric parameters have been measured, are cat¬ 
aloged in the PhotoObjAll table. This table has several 
sub-tables, known as views in SDSS: PhotoObj, which in¬ 
cludes all primary and secondary objects; PhotoPrimary, 
which includes all the primary detections or ones classi¬ 
fied as the best version of the objects; PhotoSecondary, 
which includes the duplicate detection(s); and PhotoFam- 
ILY, which includes the original undeblended source as well 
as the sources for which the deblending failed. Based on 
the morphology, the PhotoPrimary is further divided into 
Star (point sources). Galaxy (extended sources). Sky (sky 
samples), and UNKNOWN (unclassified sources). A PHOTO¬ 
PRIMARY object is classified as a point source and included 
in Star when a PSE fit provides a good approximation to its 
light profile; otherwise, it is classified as an extended sourc^ 

As the table with all the unique point source objects. Star 
is the reposit ory used by all stellar studies that use the SDSS 
Covey et al.[[2007 Juric et ar][2008[ Bochanski 


surve ' 
et al 


ey (e.g ., 
l2^ i 


including this study. Therefore, it is troublesome 
that we detected no resolved companions within 1" in our 
search. Given the scope of Star, it is of utmost importance 
to understand its completeness level. In particular, it is nec¬ 
essary to quantify the sources that were detected in SDSS but 
are missing from Star, i.e., the sources that were misclassi- 
fled. This would allow any individual study to apply relevant 
corrections for the sources that were truly missing from SDSS 
survey (e.g., faintness limit, unresolved binarity). 

To test the completeness of Star, a sample of bona fide 
stellar sources needs to be used. Fortunately, there exist 
large spectroscopic samples of confirmed stellar sources in the 

^ http://www.sdss3.org/dr9/algorithms/ 

* More specifically, a point source has PSFMag — CMODELM AG <0.145, 
where PSFMag and CMODEL are the magnitudes measured by fitting a PSE 
model and a linear combination of de Vaucouleurs and exponential models, 
respectively, for an object’s light profile. 
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SDSS catalog. We chose the DR7 catalog of 70,841 M dwarfs 
from the SDSS|West et al.|(|2011|). Each spectrum in this cat¬ 
alog was inspected by eye and verified to be a bona fide M 
dwarf with relatively high signal-to-noise. Therefore, look¬ 
ing at how these bona fide M dwarfs are cataloged in Star 
using only photometric information and how they are flagged 
should help us understand why we were unable to find bina¬ 
ries tighter than 1". As the selection algorithm for M dwarfs 
in the SDSS spectroscopic survey was not based off Star 
( |West et al.|201 Ij l, this is an ideal sample to investigate how 
many of the 70,841 bona fide M dwarfs are not included in 
Star. 


3.1.1. What is the completeness ofDR7 M dwarf sample? 

We performed a PLATE/MJD/FiBERlD-based search for the 
70,841 M dwarfs in DR7 STAR using the SDSS CasJobs por¬ 
tal and recovered only 66,001 (~93.2%). When the same 
search was performed on DR7 PhotoObj, every single spec¬ 
troscopic target had an counterpart]^ Appropriate flag cuts 
(see Section were used in both searches; when excellent 
photometry was not present, we assumed the sources were not 
real. That every single M dwarf was in PhotoObj but ^7% 
were missing in Star made it evident that they had been mis- 
classified. 

The first type of misclassification (544 instances; ~0.77% 
of all sources) happened when an M dwarf was imaged two or 
more times, and the Primary detection failed the photomet¬ 
ric flag cuts but one or more SECONDARY detection(s) passed 
them. About one-third of the SDSS footprint was imaged 
multiple times, mostly when plates overlapped each other. 
In such cases the PRIMARY / SECONDARY classification is 
based on whether the held was designated as the “primary” 
held for that area of the sky. As the quality of photometry 
plays no part in this classification, this is understandable. In 
fact, it is remarkable that only 0.77% of sources are lost; this 
incompleteness is largely negligible. 

The second kind of misclassification, which happened in 
4,295 instances, happened when the M dwarfs were classified 
as extended sources. This was either because of an nearby 
extended source contaminated and extended the light profile 
of the M dwarf or there were two partially resolved stellar 
sources. (The latter are exactly the kind of sources we are 
looking for in our binary search.) That ~6% of all stars could 
be not cataloged in Star is rather alarming. We note that this 
effect is likely exaggerated for the spectroscopic sample (as 
compared to the entire photometric sample) as a large num¬ 
ber of M dwarf spectra were acquired as they were targeted a s 
potential quasars and LRGs (|Adelman-McCarthy et al.|2006|l. 
Thus, they are more likely to be flagged as being extended 
sources. However, an incompleteness of ~6.1% is signifi¬ 
cant and should be properly accounted for studies that use the 
SDSS photometric sample. 

3.2. Photometric Distances 

Main-sequence Dwarfs: We determined distances to main- 
sequence dwarfs using photometric parallax relations that 
are calibrated from directly-measured trigonometric paral¬ 
laxes. Even though there is no spectroscopic confirmation, 
MS dwarfs dominate the stellar sources in the SDSS catalog 

^ This search was done in DR7, from which the spectroscopic catalog 
was compiled. As each DR is completely new reduction, it is expected that 
not 100% would be recovered between different DRs. Indeed, we recovered 
fewer objects when we performed the same search in DRS and DR9 catalogs. 



r-z 

Figure 2. The M subdwarfs, extreme subdwarfs, and ultra subdwarfs (purple, 
green, and red circles, respectively) are clearly separated from the M dwarfs 
(black dots). The separation is even cleaner when the photometr y is precise; 
the bottom panel shows all the M dwarfs in SDSS DR7 catalog jWest et al.| 
|201 1) while the middle and top panels only show ones with uncertainties 
<0.05 and 0.02 mags, respectively. The blue box shows the subdwaif locus 
identified in|Bochanski et al. 1(2013) for which photometric distances can be 
calculated. 


} Covey et al.|2007 1 . Three separate photometric parallax re- 
ations were used, such that the chosen magnitude and colors 
trace the stellar temp erature monotonic ally in that regime: Mg 
vs. g — i for E0-K4 (|Covey et al.|2007[ ), M^ vs. r — z for K5- 
M9, dBochanski et al.||2010[ ), and m/ vs. i — z for L0-L9 
( Schmidt et al.^2010| l dwarfs. The photometric parallax rela- 
tions are given in Table 

Magnitudes were corrected for extinct ion, with values from 
, Einkbeiner, & Da\^ 1998| l as tabulated in the 
SS database, were used in all cases. We adopted an er¬ 
ror of 0.3 mag in the calculated absolute magnitudes, which 
implies a Icr error of 14% in distance. This scatter is caused 
by a combinati on of metallicity, magnetic activity, and unre- 
solved binarity (jWest et al.|2005| |Sesar et al.|2008| |Bochan-| 
ski et al.|2011[). The presence of magnetic activity and higher 


metallicity increases the the intrinsic brightness of a star by 
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Table 1 

Photometric parallax relations used in this paper 


type 

locus 


Photometric parallax relation 

References 


F0-K4 

-0.01 < (r-z) < 0.50 

Mg 

= 2.845 -1- 1.656 (g-i) -1- 3.863 (g-i)^ - 1.795 (g-i)^ 

Covev et al. l2007l 


K5-M9 

0.50 < (r-z) < 4.53 

M,. 

= 5.190 -1- 2.474 (r-z) + 0.4340 (r-z)^ - 0.08600 (r-z)^ 

Bochanski et al. |2U1U 

L()-L9 

1.70 < (i-z) < 3.20 

Mi 

= -23.27 + 38.40 (i-z) - 11.11 (i-z)^ -|- 1.064 (i-z)3 

JSchmidt et al. i iUlU 


M subdwarf 

(g-r)> 1.50; 0.8 < r-z <2.5 

Mr 

= 7.9547 -1- 1.8102 (r-z) - 0.17347 (r-z)^ -|- 7.7038 5g-r - 1.4170 (r-z) 6g-r 

Bochanski et a . j2U 

3 

white dwarf 

|Girven et aL|J201 If 

iterative fits to DA cooling models 

Mams et ai.|jzuuo| 



as much as a magnitude |Bochanski et al.| (|2011)1. An unre¬ 
solved binary companion increases the apparent magnitude 
by as much as 0.75 mag. Based on high-resolution imaging 
studies, components of wide binaries are highly likely to har¬ 
bor close companions ( [Law et ar]|2010| l. Therefore, we are 
selecting against triple systems that contain an unresolved bi¬ 
nary. However, none of these parameters can be measured or 
even estimated from the photometry alone. Ivezic et al. ( |2008 1 
used a photometric metallicity in their photometric parallax 
relation, but that is applicable only to FGK dwarfs. No such 
relation is known for M dwarfs, except when they are sig¬ 
nificantly metal-poor subdwarfs (see below). Therefore, we 
chose to ignore the metallicity-dependence. 

M Subdwarfs: M subdwarfs are low-metallicity, main- 
sequence dwarfs typically associated with the thick disk or 
the halo. In the M spectral type, the subdwarfs can be dif¬ 
ferentiated from t he dwarfs b y th e deplet ed TiO feature in the 
optical spectrum ( |Gizis|1997| l. In |Paper l] we used the reduced 
proper motions to select subdwarf candidates and identify 70 
CPM binaries as subdwarf systems. As a photometric paral¬ 
lax relation had not been calibrated for M subdwarfs, we used 
the one for M dwarfs. As a result the distances were overesti¬ 
mated andwewerea^ identi fy only a subset of the subdwarf 
pairs. [Bochanski et al.| ( |2013| ) have since developed a statis- 
tical parallax relation using the SPSS DR7 subdwarf catalog 
( jSavcheva, West, & Bochanski |2014| l based on their redder 
g — r colors. Likely caused by increased hydride absorption, 
the redder g — r color causes the M subdwarf locus to clearly 
separate from the M dwarf locus in the (r — z, q — r) space 
( |West et al.|2004||20TT]|Lepine & Scholz|2008| ). 

Distinguishing the subdwarfs from the dwarfs, however, is 
not trivial. As shown in Figure!^ while the loci clearly sepa¬ 
rate in the {r — z, g — r) space, the tail of the M dwarf distribu¬ 
tion scatters into the subdwarf locus. Given their vastly larger 
number in the Solar neighborhood, the dwarfs overwhelm the 
subdwarf population. However, we found that the scatter is in 
large part due the uncertainty in the magnitude measure ments. 
In Figure |2| the spectroscopically confirmed M dwarfs (|West 
et al. 120111 are plotted as black dots while the spectroscopi- 
cally conhimed subdwarfs are plotted as purple (subdwarfs), 
green (extreme subdwarfs), and red (ultra subdwarfs). In the 
bottom panel, where all the st ars plotted, the dwarfs scatter 
into subdwarf locus (blue box; |Bochanski et al.|2013| ) signif¬ 
icantly. However, when only the dwarfs with uncertainties 
<0.05 (middle) and <0.02 mag (top) in the griz bands are 
plotted, the scatter decreases. While the photometric accuracy 
for SDSS is ~1% in the griz bands, imposing a 1% or 2% cut 
seems to exclude a large number of sources. This also biases 
the sample against fainter and more distant stars like the M 
subdwarfs. However, when available, exquisite photometry, 
can be used to select subdwarfs without the need for spectra. 
This could possibly be extremely beneficial to future photo¬ 


metric surveys like the Large Synoptic Survey Telescope. 

We chose to accept photometry with uncertainties <0.05 
mag to identify subdwarf candidates. We recognized that this 
will include some M dwarf interlopers and increase the rate of 
false positives. However, given the low number of subdwarf 
binaries known, this is an acceptable risk. We did not include 
the identified subdwarf systems in the statistical analysis for 
this paper. We did not include the subdwarf binaries in the 
analysis for this paper, as they are less likely to be bona fide 
systems compared to the rest of our sample. 

We calculated photometric distances to subdwarf candi¬ 
dates that fit the following criteria; 


g — r > \.S 
0.8 < r-z < 2.5 
psfMagEiTg^i^ <= 0.05 


( 1 ) 


from the Bochanski et al.| ( (2013 l relations. The uncertainty 
in the absolute magnitude is ^0.41 mag, which translates to 
a 20% uncertainty in the distance. There were over six mil¬ 
lion subdwarf candidates, as defined by Eq. in the SDSS 
photometric catalog. 


White Dwarfs: In [Paper I we used proper motions to cal¬ 
culate the reduced proper motions and identify potential WD 
candidates and search for WD companions to the low-mass 
star sample. Based on ugriz photometry alone, there is no 
conclusive way to identify potential WD companions around 
our low-mass dwarf sample. However, in the color-color 
diagrams, WDs segregate from the main-sequence stars and 


the quasars. In particular, [Girven et al. (2011|l have de¬ 


fined a (g — r, u — g) locus for hydrogen-atmosphere WDs 
(DAs) based on a sampl e of spectroscopically identified DAs 
( jEisenstein et al. 2006| l in SDSS DR7. The efficiency of 
the photometric selection for the spectroscopic sample was 
only ^62.3%, after non-DA WDs including WD-I-MS pairs 
(^8.9%), early-type MS stars and subdwarfs (~11.3%), and 
quasars (^17.2%) were removed. However, as quasars were 
specifically targeted for the SDSS spectroscopic sample, they 
are disproportionately represented. In the larger photometric 
sam ple, the contaminati on due to quasars is lik ely to be mi¬ 
nor ([Girven et al.|2011| ). Therefore, we use the [Girven et al.[ 
(|201 l[l (g — r, u — g) locus to identify potential WDs but note 
that >20% of the WD candidates will be interlopers and will 
require spectroscopic confirmation. 

We calculated the photometric distances to candidat e DAs 
by fitting the ugr iz photometry t o WD cooling mod els ( [Berg-[ 
[eron et al.|1995[ ), as specified in [Harris et al.[ ( 200^ . This al- 
gorithm hts the photometry to the model in an iterative man¬ 
ner to derive a bolometric luminosity and a distance. How¬ 
ever, the composition and mass/gravity is degenerate when 
only the photometry is available and cannot be determined. 
So we used the models with pure hydrogen atmospheres and a 
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gravity of log g = 8.0. As our adopted (g — r,u — g) locus se¬ 
lected only DAs, the composition should introduce significant 
uncertainties in the distances. However, incorrect distances 
will be derived for WDs with unusually high mass/gravity 
(^15% of all WDs), unusually low mass/g ravity (^10% of al l 
WDs), or helium-dominated atmospheres (|Harris et al.|2006|l. 

A comparison of our photometric distances and the spectro¬ 
scopic distances for th e DA WDs in the SDSS DR7 catalog 
( |Kleinman et al.|2013| l showed a 14-20% scatter (J. Andrews, 
private communication). To be conservative in our matching 
process, we adopted 14% as the error in our photometric dis- 
tan ces. There we r e 56,5 05 WD candidates, as selected using 
the IGirven et al.| ( |201 1| | (g — r, u — g) locus, in the SDSS 
photometric catalog. 


3.3. Binary Candidate Selection 

We searched for stellar sources that had been classified as 
Primary detections within angular separations of 1-20" of 
our sample of 24,036,982 low-mass stars. The inner search ra¬ 
dius of 1" was determined by the SDSS database, as discussed 
in Section 3.1 We chose the outer search radius to be 20", af¬ 


ter which the number of chance a lignmen ts for pairs without 
proper motions were significant in Paper I| (see Figure[T]). The 
search was conducted using in the NEIGHBORS table in the 
SDSS CasJobs Query interface. Photometric quality cuts, as 
described above, were performed. We then matched the pho¬ 
tometric distances id) to within la. However, as the error in 
the distances is a percentage error, we also required the dif¬ 
ference in the distance to be less than 100 pc to be classified 
as a candidate binary pair. Thus, our binary candidate pairs 
matched the following criteria: 

0 = 1 - 20 " 


A d < min(lcrA d, 100 pc) 

d < 2500 pc (2) 

| 6 | > 20 °. 


All stars were selected to have good photometry and PSF- 
MagErr < 0.10 mag for the bands that are used in their 
selection and analysis. We identified 514,424 dM-|-MS, 1212 
sdM-fsdM, and 642 WD-PdM candidate pairs. 


3.4. The Galactic Model: Assessing False Positives in the 
Binary Candidate Sample 


Despite the rigorous nature of our candidate selection, chance 
alignments will be present in a sample of wide visual bina¬ 
ries. Such chance alignments arise from the measurement un¬ 
certainties in the parameters used in the selection criteria, as 
well as via the th e inherent spreads in these parameters in the 
Galaxy ( [Paper I| l. The number of chance alignments grows as 
a function of the angular separation (cx 9^) and distance. As 
our search does not include any kinematic matching, the prob¬ 
ability of chance alignment is particularly high and, there¬ 
fore, needs to be rigorously assessed for each and every can¬ 
didate binary pair. Such a quantit ative ass essment sifts out 
false positives from the sample. In [Paper I| we built a Monte 
Carlo-based Galactic model that recreated the stellar popu¬ 
lations along the line-of-sight (LOS) of a candidate binary 
and calculated the probability of chance alignments. Based 
on empirically-measured parameters for the Milky Way, the 
model accounted for the variations in the stellar number den¬ 
sity and space velocities, which become important beyond the 


Table 2 

Galactic Structure Parameters 


Component 

Parameter name 

Parameter description 

Adopted Value 


P (-^0,0) 

stellar density 

0.0064 


/thin 

traction** 

i "/thick "/halo 

thin disk 


scale height 

260 pc 



scale length 

2500 pc 


/thick 

fraction^ 

9% 

thick disk 

-^thick 

scale height 

900 pc 


-^thick 

scale length 

3500 pc 



fraction" 

0.25% 

halo 

^halo 

density gradient 

2.77 


q (= c/a)'° 

flattening parameter 

0.64 


Note. — The parameters were measured using M dwarfs for the disk jBochanski| 
|et al.|20fQ) and main-sequence turn-off stars for the halo (June et al.|2Q08) in the SDSS 
footprint. 


^ Evaluated in the solar neighborhood 
^ Assuming a bi-axial ellipsoid with axes a and c 


Solar Neighborhood. We employed the same model to assess 
the probability of chance alignment for the binary ca ndidates 
identified here. The model was described in detail in [Paper I{ 
here we only provide a brief synopsis. 

Instead of the computationally implausible task of simu¬ 
lating the entire Galaxy with > 10^^ stars, we recreated a 
30" X 30" cone centered at the (a, 5) of each primary out to 
a distance of 2500 pc. First, we calculated the total number 
of stars in the conical volume by integratin g stellar number 
densit y profiles that assume a bimodal disk ([Bochanski et al.[ 
2010[l and an ellipsoidal halo Puric et al.|2(j08[l. ITie model 
parameters are given in TablelS] While a two-component 
model is an oversimplification oTthe com plicated scale height 
distribution ( Bovy, Rix, & Hogg||2012| , it is easy to model 
and suffices for our purpose of recreating a random Galaxy. 
Each LOS was then repopula ted with stars using the rejec¬ 
tion method ([Press et al.[1992jl. The rejection method ensured 
that the stars were randomly redistributed while following the 
overlying stellar number density distribution function. Each 
star that is generated has a a, 5, and distanca^ The total num¬ 
ber of stars in an LOS ranged between 2000Mt)000, enough to 
both produce small-scale density variations as are seen in the 
Milky Way and perform calculations with relative ease. How¬ 
ever, as the number density profiles are smoothed functions 
that were made to fit to the Galactic field, our model does not 
produce larger density variations like moving groups, open 
clusters, or streams. 

With the simulated galaxy, we then counted how many stars 
were found in the ellipsoid that is centered at the a and d 
of the binary candidate and defined by its angular separation 
and its distance errors. This was the same criteria that we 
used in the search for candidate systems in the SDSS data. 
Indeed, to conduct the exact same search, we chose to not 
convolve the luminosity or mass function into the model and 
rather searched for all stars. As all stars in our simulation 
are single stars, any star that satisfied the search criteria is a 
chance alignment. Therefore, the average number of chance 
alignments is an assessment of the probability that our candi¬ 
date binary is a false positive. For each candidate pair, we ran 
1000 Monte Carlo realizations. If the probability of chance 
alignment, Pf < 0.10, we ran a further 4000 realizations for 
better resolution. 


InjPaper ijwe modeled the 3D velocities for each star, but we do not do 
so here as we are not utilizing proper motions 
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Figure 3. A gri composite collage of SLoWPoKES-II binaries, 48” on a side. We have identified 105,537 bona fide binaries, with at least one 
component later than K5, without using proper motions. 


4. RESULTS: THE SLoWPoKES-II BINARY SAMPLE 

Following |Paper I| we classified candidate pairs with a prob¬ 
ability of chance alignment, Pf < 0.05 as real binaries. We 
note that this limit does not have any physical motivation 
but was chosen to minimize the number of spurious pairs. 
This cut results in 105,537 dM+MS, 450 WD+dM, and 944 
sdM+sdM binary systems with separations of 1-20". 141 of 
dM+MS binaries are VLM binary candidates Table 4, with 
i — z colors redder than the median M7 dwarf for both com¬ 
ponents. Despite the lack of kinematic information, we dub 
this catalog SLoWPoKES-II. The data for the parrs are tabu¬ 
lated in and is also available in an online visualization portal^. 

This represents a significant increase over the SLoWP oKES 
catalog of 1342 CPM binaries we presented in [Paper I| Each 
binary has a low probability of chance alignment, with the 
threshold set at Pf < 5%. Using that threshold, we might 
expect 5277 of the pairs to be false positives. However, each 
binary has its own probability of chance alignment, based on 
the size of the binary and its position in the Galaxy, that in 
most cases is significantly lower than 5%. Based on those 
probabilities, we expect only 2464 (or 2.33%) of the binaries 
to be falsepositives. 

Figure l^shows a collage of gri composite images, 48" on 
a side, ofl5 SLoWPoKES-II binaries. The inferred spectral 
types of MS dwarfs, based on their r — z colors, are shown 
for each component. Table 3 summarizes the SLoWPoKES- 
II catalog, with the properties of the systems and positions 
and photometry of both components included. Table 4, Ta¬ 
ble 5, and Table 6 summarize the candidate VLM, WD+dM, 
sdM+sdM systems, respectively. While these systems are 
also analyzed with the Galactic model and required to have Pf 


< 5%, their classification as an sdM, VLM, or WD is based 
on photometry alone. This adds a level of uncertainty that the 
Galactic model cannot quantify. For that reason, these sam¬ 
ples are likely to be contaminated at a higher rate than the 
SLoWPoKES-II sample in Table 3. Thus, we do not include 
samples in the analysis presented in this paper. 

Both the SLoWPoKES and SLoWPoKES-II catalogs are 
publicl y available on the w orldwide web® via the Filtergraph 
portal (Burger et al. 20131. We hope this will enable to the 
entire community to interact with and visualize this large data 
set using a dynamic, multi-dimensional plotting applet. The 
table view Filtergraph is also an easy medium to select targets 
for followup observations. 


4.1. Identifying bona fide wide binaries without kinematics 

Figure shows the number de nsity of the 514,424 binary 
candidates, identified in Section |3.3| as a function of the 
probability of chance alignment, Pf, and angular separation. 
As would be expected, the number of candidate pairs grows 
sharply with angular separation. There is also a large spread in 
Pf at a given separation, reflecting the varying stellar densities 
along different lines-of-sight. Especially at the large angular 
separations, the Pf is quite high for a large number of candi¬ 
dates, indicating that they are chance alignments. Therefore, 
we applied the Pf < 5% cut, shown as the red line, for the can¬ 
didates to be included in the SLoWPoKES-II catalog. This is 
a fairly conservative cut but serves to minimize the number of 
false po sitives. It is also the same threshold that we used in 
Paper I| although the SLoWPoKES systems had proper mo¬ 
tions and, hence, should have a much lower rate of false pos¬ 
itives. The middle panel shows the histogram of fraction of 
candidate binaries that passed this threshold as a function of 
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Figure 4. The number density, as represented by the contours, of the 514,424 
binary candidates as a function of Pf as calculated by the model vs. the angu¬ 
lar separation and distance. The red line is the threshold Pf that we adopted 
for the SLoWPoKES-II catalog. The fraction of candidates that passed the 
Pf threshold—which quantifies the number of chance alignments expected 
within the bounds of that particulai* binary—are shown in the middle panels 
while the distribution of the resultant 105,537 SLoWPoKES-II binaries are 
shown in the bottom panels. 

angular separation. All of the candidate pairs within 3" but 
only ~2% of the pairs at 20" were vetted to be real binaries by 
the Galactic model. The bottom panel shows the angular sep¬ 
aration histogram for the 105,537 binaries in SLoWPoKES-II 
catalog. This distribution is strongly skewed towards smaller 
separations, with 51% of the pairs at l"-7" and 72% at 1"- 
10 ". 



Figure 5. The fraction of binaiy candidates that meet the Pf < 1%, 5%, 10%, 
and 50% thresholds as a function of angular separation. The number of pairs 
that meet each threshold is also shown. We have chosen Pf < 5% (dark line) 
for the SLoWPoKES-II catalog, resulting in 105,537 wide binaries. 


Figure|^shows the distribution of the fraction of binary can¬ 


Figure 6. The number density, as represented by the count-ours, of the 
514,424 binary candidates as a function of Pf and distance. The red line is the 
threshold Pf that we adopted for the SLoWPoKES-II catalog. The fraction of 
candidates that passed the Pf threshold are shown in the middle panels while 
the distribution of the resultant 105,537 SLoWPoKES-II binaries are shown 
in the bottom panels. 


didates that pass various Pf thresholds; < 1%, 5%, 10%, and 
50%, with the total numb er of can didates also shown. Obvi¬ 
ously, as we discussed in |Paper I| the choice of a threshold 
is subjective and depends on the need. For example, if a fol¬ 
lowup study needed to minimize the number of false positives, 
a sample with Pf < 1% should be chosen. However, if the 
study needed bright, rare binaries (e.g., WD-fdM binaries), it 
might need to tolerate the high number of false positives and 
select candidates with Pf < 10%. With such a large number 
of binary candidates, there is a flexibility in the number and 
types of binary systems that is available to followup studies. 

Figure]^ shows the number density of all candidate pairs as 
a function of Pf and photometric distance. In general, Pf is 
flat as a function of the distance, indicating that the distance 
of a given pair has no significant effect on how likely it is 
to be a chance alignment. More significantly, the number of 
candidate pairs peaks around ^800 pc and declines smoothly 
afterwards. This is because an additional selection criteria, 
A d < 100 pc (Eq.|2]), which becomes effective at large dis¬ 
tances and rejects pairs with A d within 1 ciAd but larger than 
100 pc. So even as the search volume grows larger, the num¬ 
ber of candidates actually decreases. However, the resulting 
candidates are more likely to be binaries, as shown in the mid¬ 
dle panel; the fraction of candidate binaries that pass the Pf 
< 5% threshold reaches a minimum at ~700 pc and increases 
sharply afterwards. In fact, the acceptance rate for binaries at 
^2200-2500 pc is higher than for binaries at any other dis¬ 
tance. At these distances, A d < 100 pc implies that the dis¬ 
tances match within 6-7% of each other, explaining the low 
rate of false positives. Therefore, as shown in the lower panel 
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of Figure ^ there are a relatively large number of identified 
binaries atlarge distances. 

4.2. Distribution of binary separations 



Figure 7. The distribution of projected physical separations for the 
SLoWPoKES-II binaries (top) and the SLoWPoKES CPM binaries 
(bottom). As SLoWPoKES-II was restricted to 9 < 20", it only 
probes binaries smaller than ~50,000 AU, thus missing the widest 
systems. Even so, the slowpokes-II distribution cannot be fit by a 
single functional form. We fit two polynomials with a break at a = 
10^'AU, shown in dashed red and solid blue lines in the figure. 
This inflection point suggests the presence of multiple populations 
of wide binaries. 

Figure shows the distribution of projected physical sep- 
arations of the SLoWPoKFS and SLoWPoKFS-II binaries IH 
The SLoWPoKFS distribution exhibited a bimodal distribu¬ 
tion, which we interpreted as either the presence of two for¬ 
mation modes or the preferential d estructio n of the widest 
pairs in the timescale of a 1-2 Gyr (|Paper I||. These binaries 
were 7-180" systems, within an average distance of 1000 pc. 
In SLoWPoKES-II, we are not sensitive to binaries as wide, 
as the lack of proper motions limits us to angular separa¬ 
tions of 1-20"; however, the binaries are up to 2500 pc away. 
Clearly, the selection biases and the resultant distributions of 
the identified binaries involved are markedly different. For 
example, the steep falloff of the SLoWPoKES-II distribution 
at a > 10"^'^ AU is due to our insensitivity to binaries at those 
separations. 

*' In|Paper l|we used a statistical conection to scale the projected physical 
separation to semi-major axis; however, |Dupuy & Liu| (2011) have shown 
that the correction factor is significantly smaller and highly dependent on the 
orbital parameters. They calculated correction factors based for their sample 
of <4 AU binaries. The sample presented here are in a completely different 
separation regime, making it unlikely that the same correction factors can be 
used. Therefore, we have eschewed semi-major axis and chosen to work in 
projected physical separation space. 


In Paper I we found a dip in the distribution of binary sepa¬ 
rations at a ~ 10"^ AU (bottom panel of Figure^, but at first 
glance this is not apparent in the new SLoWPoKES-II distri¬ 
bution (top panel). However, there is an inflection point at 
a ~ 10^'^ AU (2500 AU), where the slope of the distribution 
becomes noticeably shallower. In Figure two power-law 
fits are shown on either side of the inflection at a ^ 10^'^ AU. 
This inflection is at a similar separation scale as seen in pre- 



common, as was originally predicted by Opik (19241. How¬ 
ever, it is much smaller than the bimodahty seen in Paper I| at 
^20,000 AU. Whether present at formation or sculpted later 
in life, the multiple modes in binary populations is intriguing 
and suggests that interstellar interactions are more common 
than thought. 

4.3. Binary mass distribution 

Figure shows the g — i and r — z color distributions for 
the primary (left) and secondary (right) components of the 
SLoWPoKES-II binaries. The primary is defined as the com¬ 
ponent with the bluer r — z color. The histograms for the 
color distributions are shown along the top and sides of each 
panel. The inferred spe ctral types, based on the median color- 
spectral type relations ( |Covey et al. 120071 West et aTl20111) , 
are shown along the top axes. In general, the SLoWPorCFS-II 
binaries reflect the SDSS low-mass dwarf population in both 
g—i and r — z colors. The distributions peak between the M2- 
M4 spectral types, same as the held mass function (|Bochan- 
ski et al.|2010|l. The primary color distributions also exhibit 


a pileup at the bright end of the distribution, which is most 
likely a selection bias. SLoWPoKES-II binaries span a large 
range of colors (and masses). There are a significant num¬ 
ber of primary components from early-mid G to the mid-M 
dwarfs while the secondaries extend from K5 to the M9 spec¬ 
tral types. The blue end of the secondary component distri¬ 
butions are defined by the color cuts imposed on our initial 
target sample (r — i >0.3, i — z > 0.2). While SLoWPoKES 
contained only a handful of binaries after ~M6-M7, there are 
now a significant number of binaries at the low-mass end of 
the stellar main sequence at large distances. 

Figure 1^ shows mass ratio distribution of SLoWPoKES-II 
binaries as a histogram (top) and a cumulative distribution 
function (bottom). The black line and dots represent the en¬ 
tire catalog while the red, purple, blue, and cyan show distri¬ 
butions for different primary masses. The masses were cal¬ 
culated by inte rpolating from their g — i, r — z colors, or 
i — z based on Kraus & Hillenbrand] ( |2007a| l. We set the 
lowest mass at 0.075 Mp, the hydrogen burning minimum 
mass ( [Burrows et al.||T997[ ), as the i — z-mass relation is 
even less constrained at lower masses. All of the distributions 
are skewed towards similar masses, with notable deficit at the 
lowest mass ratios, below q ~0.4. This is mostly due to an ad¬ 
ditional selection criterion that at least one component be K5 
or late (<0.70 M©), which means the lowest mass ratio possi¬ 
ble for a K5 primary is ^0.11. Moreover, the steep mass- 
luminosity ratio among the M spectral types means SDSS 
photometry would only detect a M6 or earlier companion for 
the K5 primary, restricting the lowest possible mass ratio to 
^0.17. Similarly, for a M4 primary (~0.20 M©), the ob¬ 
served mass ratio is always greater than 0.39. After q ~0.4, 
the various cumulative distributions show a relatively smooth 
progression that is similar for the different masses. There is a 
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Figure 8. The g — i and r — z color distributions for the primary (left) and secondary (right) components of the SLoWPoKES-II binaries, with the contour levels 
at 10, 100, 1000, and 5000. The histograms of the distributions are plotted along the sides, and the inferred spectral types are also shown. We defined the primary 
component to be have the bluer r — z color. Both the primary and the secondary distributions show definitive peaks at ~M2-M4 spectral types in both the g — i 
and r — z colors. As our initial target sample used a cutoff on the r — i and i — z colors at ~K5 spectral type, there is a sharp cutoff at the blue end of the r — z 
distribution for the secondary component. 


preference for a larger q at the lowest masses and a smaller q 
at the highest masses, but those are most likely due to selec¬ 
tion biases as discussed above. 

5. DISCUSSION 

We have assembled a large sample of 105,537 wide binaries 
with projected physical separations of ~ 1000-60,000 AU. 
The identification was based on matching 3D position of stars 
in the SDSS photometric sample and assessing the proba¬ 
bility of chance alignment with our Galactic model. All of 
the binaries in the SLoWPoKES-II catalog have a probabil¬ 
ity of chance alignment of <5%. For comparison, the origi¬ 
nal SLoWPoKES contained 1342 CPM binaries with separa¬ 
tions of 7"-180" (500-100,000 AU) at -100-800 pc. While 
our initial aim was to combine these two samples for a de¬ 
tailed analysis, the regimes they probed are much differ¬ 
ent. For example, despite extending to >2000 pc, the widest 
SLoWPoKES-II binary is a factor of two smaller than the 
widest SLoWPoKES binary. Proper motions allowed for the 
identification of CPM binaries up to 180", whereas we were 
restricted to 20" in this paper. Thus, even though the two sam¬ 
ples are are neither complete nor continuous, they are comple¬ 
mentary and are both of high fidelity. 

5.1. SLoWPoKES-II: A rich, diverse catalog of binaries 

In addition to being the largest catalog of wide binaries, 
SLoWPoKES-II contains a diversity of systems—in mass, 
mass ratios, metallicity, separations, evolutionary states, and 
Galactic positions—that should facilitate followup studies to 
characterize the properties of low mass stars. In addition, we 
have identified three specific subsets of candidate binaries that 
are rare: (i) wide VLM systems, (ii) WD+dM systems, and 
(iii) sdM-fsdM binaries. They could prove uniquely useful in 
characterizing those stellar types. 

The diversity of systems should allow for in-depth probes of 
different aspects of low-mass stellar physics. Specifically, the 
flexibility offered by the large sample size and all-sky nature 


of SLoWPoKES-II could be exploited by large spectroscopic 
surveys where fibers are often available in certain parts of the 
sky. Eor example, we were awarded 1000 fibers in the B aryon 
Oscillation Spectroscopic Survey of SDSS-III (BOSS; |Daw-| 
son et al.l2013|) to acquire spectra of 500 SLoWPoKES-II bi¬ 
naries. Inis was in addition to the 1000 fibers awarded ear¬ 
lier for the SLoWPoKES binaries. We are using the com¬ 
bined sample —1000 binaries to ascertain the age-activity and 
metallicity relationships in low-mass M dwarfs (Massey et al. 
in prep.). 

5.2. Multiple pathways for wide binary formation 


In Paper I we noted the presence of a bimodality in the 
physical separation distribution of the binaries (Eigure |7]). 
When compared to t he dynamical dissolution timescales 
( Weinberg et al.|[l987 1 , the bimodality in the SLoWPoKES 
distribution led us to suggest that it comprised of two differ¬ 
ent population of binaries with a break at —20,000 AU: (1) 
a “wide” population that is dynamically stable over —10 Gyr 
and (2) an “ultra-wide” population of young, loosely-bound 
systems that will dissipate in 1-2 Gyr. In a followup high- 
resolution imaging study with the Keck II and Palomar Laser 
Guide Systems with Adaptive Optics, we found that the fre¬ 
quency of a cl ose companion wa s higher in wide binaries than 
in single stars ( Law et al.|2010| l. Moreover, the frequency in¬ 
creased significantly with wide binary separation, suggesting 
that different formation modes were at work. 

In SLoWPoKES-II we have found further evidence of mul¬ 
tiple populations in the distribution of physical separations 
(Eigure |7]i. While more subtle than in |Paper I| a break is 
clearly evident at — 2500 AU . This is much smaller than the 
critical separation in Paper I| but approximately where we saw 
the first peak in the SLoWPoKES sample. It is also at the 


same separation scale as seen in previous studies (Opik 1924 
Allen et al.||2000[ |Lepine & Bongiorno||2007|l. That an in 
flection in the sep^ation distribution stands out amongst the 
myriad of selection biases and incompletenesses is quite sig- 
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nificant. Further investigation into these multiple populations 
and their origins is clearly warranted. Whether wide binary 
distrib utions are segregated at birth (e.g., |Reipurth & Mikkola 
2012|l or sculpted as they traverse around the M ilky Way (e.g., 
Weinberg et al.|19^ Jiang & Tremaine|2010 l has important 
implications on our interpretation of observed binary popula¬ 
tions and, consequently, on our understanding of binary star 
formation. 

There are no signatures of multiple populations in the color 
(a proxy for mass; Figureor mass ratio (Figure]^ distri¬ 
butions. Fven if there were such signatures were present, 
identifying them from amongst the selection biases without 
further data would be extremely difficult. Recently, there has 
been extensive discussion on the formation and stability of the 
extremely wide systems. This has partly been motivated by 
the large samples of wide binaries that have become available 
over the past decade (e.g., Chaname & Gould|2004[|Lepine & 


Bongiorno|2007 [ Sesar et al.|2008[ Dhital et al.|2010| l. In ad¬ 

dition, such wide binaries have been identified at very young 


_ very young 

ages in Orion (|Connelley, Reipurth, & Tokunaga|2009|l and in 
Taurus and Upper Sco ( |Kraus & Hillenbrand|2009b| l. Numer- 
ical simulations have explored mechanism that would form 
wide binaries primordially or via dynamical interactions, with 
three different pathways suggested. First, observations of pro- 
tostars at the end of long, extended filaments of molecular gas 
(Tobin et al.||2010]l argues that wide binaries can form with 



0.0 0.2 0.4 0.6 0.8 1.0 

Mass ratio (q = M^/Mj) 

Figure 9. The mass ratio distribution of SLoWPoKES-II systems shown as 
histogram (top) and a cumulative distribution function (bottom). The black 
line and dots show the all the SLoWPoKES-II binaries while the red, pur¬ 
ple, blue, and cyan show the distributions for systems with different primary 
masses, mas ses were calculated by i nterpol ating from their g — i and r — z 
colors using |Kraus & HillenbrandH2007af; the minimum mass was set at 
0.075 M0. There is a large bias towards systems with similar masses, with 
the median mass ratio being ~0.78. While the systems with larger primary 
masses seem more likely to have low mass ratios as compared to systems 
with low-mass primary, it is largely due to biases inherent in SDSS and our 
identification techniques. 


primordial separations of ~0.1 pc. Second, primordial triple 
or quadruple systems become extremely hierarch ical by scat¬ 
tering or even ejecting one of the components ( [Reipurth &| 
|Mikkola|2012] l. Third, simulations have shown wide binaries 
forming in the halo of clus ters via dynamical interactions or 
as they escape the cluster (Bate & Bonnell 2005 Kouwen- 


hoven et al.|2010HMoeckel & Bate|2010HMoeckei & Clarke 


201 l|l. These pathways are unique and expected to imprint 
distinct signatures on the resultan t population. For exam¬ 
ple, the jReipurth & Mikkola \20\2\ pathway would produce a 
high frequency of triples and quadruples while dynamical in¬ 
teractions would produce a very high proportion of low-mass 
binaries. 

There is no clear evidence from the SLoWPoKFS samples 
that any of these pathways are dominant. Rather, there is rea¬ 
son to believe all of the pathways could be operational. With¬ 
out a reliable age indicator there is no way to observationally 
determine if components of an wide binary formed together. 
In fact, it might be impossible to do so even with an age in¬ 
dicator as the age difference of a few million years can be 
indistinguishable once the stars are in the Galactic field. We 
suggest that there is no principal pathway for wide binary for¬ 
mation. Determining the relative efficiencies of the different 
formation modes would be more valuable. 


5.3. Wide VLM binaries: Too numerous for the ejection 
hypothesis ? 

SLoWPoKF S did not contain wide very low-mass (VLM; 
M < 0.1 Mq Burgasser 20071 or brown dwarf (BD) bina¬ 
ries. They were too red and too faint to be detected in the 
USNO-B survey and, therefore, did not have measured proper 
motions. Finding more VLM/BD systems was one of our pri¬ 
mary motivations for identifying binaries without proper mo¬ 
tions. Only eleven wide (>100 AU) VLM/BD system s are 
currently known, i ncluding a VLM triple at 820 AU (|Bur- 
gasser et al.||2012|l and binaries at 5100 AU (jArtigau et al. 

(jRadigan et al.|2009]l. With the formation 
processes and techniques to measure properties for VLM/BDs 
not completely understood, the value of a larger sample could 
not be overstated. 

Until recently, the p rocesses by whic h very low mass 


(VLM; M < 0.1 M© Burgasser 


2007)1 stars and brown 


dwarfs (BDs) form were thought to be different from those 


I 


of higher-mass stars. Reipurth & Clarke ( 2001| ) suggested 
that VLM/BDs are protostars that are ejected from their natal 
cores before they can accrete and grow to stellar m asses. Such 
ejections were seen in numerical simulations (e.g., [Bate et al. 
2002 jl, and all observed VLM/BD binaries had enough bind- 
ing energy to have survived such ejections ([Burgasser et al.| 
2TO3| Close et al. 12003 i. However, with the discovery of each 
subsequent VLM/BD binary wider than ~100 AU, the viabil¬ 
ity o f the ejection hypothesis had been calle d into question 
(e.g., Dhital et al.|201 1[ Burgasser et al.|20T2 1 . 

Figure [10[ shows theprojected physical separation distribu¬ 
tion for a sample of 14 1 binaries with i— z > 1.14, the median 
color of an M8 dwarf ( West et al.|201 l[ l. While spectroscopic 
followup are ongoing to confirm their spectral types (a proxy 
for mass), it seems clear that significant number of systems 
with separations of 1000-5000 AU exist in the Galactic field. 
We note that photometric colors can be particularly deceiv¬ 
ing in the VLM/BD regime, so the need for spectra cannot be 
overemphasized. However, such a large number of in wide 
VLM binary candidates can inform us a lot about their forma- 
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Figure 10. Distribution of SLoWPoKES-II wide binaries, at the end of the 
main sequence. Thes e were selected by requiring i — z > 1.14, the median 
color for a M8 dwarf JWest et al.|201 1). As only eleven wide VLM binaries 
are currently known, a larger population has strong implications on how we 
interpret star formation at the lowest masses. 

tion. These systems most definitely did not form primor- 
dially via ejection, as their binding energy is too small to have 
survived a dynamical kick. With 11 previous systems and 141 
candidates presented here, we cannot explain away the wide 
VLM binary population as exceptions. 

The alternative pathway to VLM/BD formation is that they 
form in a manner similar to stars via gravoturbulent fragmen- 


tation (Hennebelle & Chabrier|201 1| Jumper & Fisher| 

b"; 


The detection of the first pre-BD core, Oph B-11 (| Andre et al. 
2012|l, has helped further that idea significantly. Our sample 


does not allow for rigorous testing of this hypothesis; we will 
wait for followup observation that will better characterize the 
binaries. However, neither the separation distribution nor the 
mass ratio distribution show any significant change as a func¬ 
tion of mass. While the systems with less massive primaries 
are more likely to be equal-mass (Figure]^. This characteris¬ 
tic is a bias of our sample mostly due to insensitivity to lower- 
mass, fainter companions and to setting the minimum mass at 
0.075 M©. The trend of more equal-mass binaries is also a 
gradual change at all masses, with no significant break at the 
VLM/BD regime. 

There is no suggestion of multiple populations of wide 
VLM binaries, as was observed for wide stellar binaries in 


Section 5.2 However, we cannot discount the possibility 
that these wide VLM binaries were bound post-formation 


( Kouwenhoyen et al.|2010[|Moeckel & Bate|2010HMoeckel &| 

|Clarke||201 l| l. Numerical simulations show that a molecular 

core typically fragments into 3-5 cores, whence the smaller 
cores are ej ected before they can grow into stellar masses 
( Bate| 20T^ . Open clusters could have hundreds of these 
VLM protostars floating around, increasing the likelihood 
of multiple ones interacting at the same time and forming 
wide, gravitationally-bound systems. The efficiency of these 
capture-like processes could be much higher in the VLM/BD 
regime. 


6 . CONCLUSION 

We have identified the SLoWPoKES-II catalog of 105,537 
wide, low-mass binaries without using proper motions. While 
false positives are inherent in all statistical samples, we have 
required a relatively stringent probability of chance align- 

For reference, the current sample of all VLM binaries is >120 (T. 
Dupuy, pnv. comm ). 


ment, Pf, as calculated by our Galactic model to be <5% for 
each of our binaries. Most binaries have a much smaller prob¬ 
ability of chance alignment. The entire sample is expected to 
have 2464 (2.3%) chance alignments. In addition, we have 
identified 944 sdM-bsdM and 450 WD-l-dM candidate bina¬ 
ries, which also have a <5% probability of chance alignment. 
However, as their identification as an sdM or WD is based on 
photometry alone, there is a higher probability of them being 
false positives. 

The STAR table in SDSS suffers from a ~6% incomplete¬ 
ness. In particular, ^ 6% of the spectros copically confirmed, 
bona fide M dwarfs (|West et al.||2011|l are missing are not 
included in STAR. Instead, they are classified as extended 
sources due to the presence of a nearby or partially resolved 
stellar object. This incompleteness could be exaggerated for 
the spectroscopic sample as a significant number o f M dwarfs 
were observed as potential galaxies and quasars ( Adelmaru] 


l.|2010| l. 
elv to ha 


This is 
have been 


McCarthy et al.||20061l. We advise caution when using the 
STAR table, especially for studies that seek to do a complete 
census or build a complete sample of stars in the Galaxy. 

SLoWPoKES-II binaries have projected physical separa¬ 
tions of 500-50,000 AU at distances up to 2500 pc. We are 
not se nsitive to binaries as wide as the ones identified in B 
per Ijl as only searched within 20". Given the dynamic mag- 
nitude limits of the SDSS, our sample is likely largely dom¬ 
inated by a confluence of selection biases making is hard to 
construct or even complete samples. For example, the color 
distributions for both the primary and secondary components 
show a higher proportion of G- and K-type dwarfs, as com¬ 
pared to the held population ( Bochanski et al.[ 
likely because they are brighter and more likely 1 
detected. With complete samples not feasible, the main con¬ 
tribution of SLoWPoKES-II is the large, diverse sample of 
binaries that will facilitate detailed studies of low-mass star 
properties. Sub-samples are already part of various studies 
to probe the age-activity, rotation-age, and metallicity of M 
dwarfs. Most importantly, SLoWPoI^S-II contains a signifi¬ 
cant number of binaries at the late-M spectral types, enabling 
the studies to probe into the very bottom of the main sequence. 

The distribution of binary physical separations exhibits a 
marked inflection which can be interpreted as repr esenting 
two binary populations, consistent with our results in|Paper I| 
and|Law et al.|P0101l. While the incompleteness in our sam¬ 
ple and lack of any way to quantify or correct for that incom¬ 
pleteness prevents us from characterizing the populations, the 
separation distribution exhibits a marked inflection which can 
be fit by two more polynomials. Combined with our results 


Paper I| and Law et ah] (2010 1 , we infer that the wide bi- 


nary population in the Galactic field is composed of multiple 
different populations, either from different formation mode or 
via different dynamical histories. 

We have identified 141 wide binaries in which both compo¬ 
nents are VLMs. This is 7x larger than the current sample of 
wide, VLM binaries. While spectroscopic data are needed to 
confirm their VLM status, it is becoming clear that wide VLM 
binaries are not exceptions. These wide systems are critical 
in understanding VLM/BD formation, as their binding ener- 
gies are too low to have formed via the ejection mechanism 
( Reipurth & Clarke|2001| l. As searches go deeper and identify 
more wide, VLM binaries, it looks unlikely that VLM/BDs 
are formed primarily via ejection. The data indicate that there 
is no change in formation mechanism with stellar mass, sug¬ 
gesting VLM/BDs are also formed via gravotur bulent frag- 
mentation, like their more massive counterparts ([Hennebelle | 
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I catalogs are avail¬ 


able on the Filtergraph portal^. (Burger et al. 2013 1 . The 
portal allows for dynamic plotting of these large data sets and 
real-time filtering of the data using user-specified criteria. It is 
useful for target selection as its smoothly transitions between 
graphical and tabular visualizations of the data and allows for 
selection of data points with the keyboard or the mouse. 

Lastly, we have demonstrated a methodology to identify 
bona fide wide binaries using just photometry together with 
a galactic model to assess fa l se pos itives. Really, this tech¬ 
nique harks back to |Michell| (|1767|l’s discovery of “double” 
stars, where he argued that some pairs of stars were too close 
to each other, as compared to mean distances between stars, 
to be unrelated. We have used the all-sky data from SDSS and 
our understanding of the Galaxy’s stellar distribution to statis¬ 
tically identify bona fide binary stars without the advantage of 
proper motions. 
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